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ABSTRACT 

We present new results on the frequency distribution of projected H I column densities f(Nm,X), 
total comoving covering fraction, and integrated mass densities phi of high redshift, H I 'disks' from a 
survey of damped Lya systems (DLAs) in the Sloan Digital Sky Survey, Data Release 5. For the full 
sample spanning z = 2.2 to 5 (738 DLAs), f(Nm, X) is well fitted by a double power-law with a break 
column density N d = io 21 - 55±0 04 cm~ 2 and low/high-end exponents a = -2.00 ± 0.05, — 6.4j^g. The 
shape of f(Nm, X) is invariant during this redshift interval and also follows the projected surface 
density distribution of present-day H I disks as inferred from 21cm observations. We conclude that 
H I gas has been distributed in a self-similar fashion for the past 12 Gyr. The normalization of 
f(Nm,X), in contrast, decreases by a factor of two during the ~ 2 Gyr interval from z = 4 to 2.2 
giving corresponding decreases in both the total covering fraction and pni- At z w 2, these quantities 
match the present-day values suggesting no evolution during the past « 10 Gyr. We argue that the 
evolution at early times is driven by 'violent' processes that removes gas from nearly half the galaxies 
at z w 3 establishing the antecedants of current early- type galaxies. The perceived constancy of phi, 
meanwhile, implies that H I gas is a necessary but insufficient pre-condition for star formation and 
that the global star-formation rate is driven by the accretion and condensation of fresh gas from the 
intergalactic medium. 

Subject headings: galaxies: evolution — intergalactic medium — quasars: absorption lines 



1. INTRODUCTION 

In the current paradigm of galaxy formation within 
CDM cosmology, baryons accrete, dissipate, and settle 
to the centers of dark matter halos. Gas with a non- 
negligible angular momentum will form an H I 'disk' 
with typical surface densities exceeding 1M Q pc -2 (or H I 
column densities, Am > 10 20 cm~ 2 ). Various processes 
(e.g. merger induced shocks, secular evolution) inspire 
the formation of molecular clouds that cool, fragment, 
and initate star formation. Finally, stellar feedback (e.g. 
winds, supernovae), AGN activity, galaxy interactions, 
and even magnetic fields (Wolfe ct al. 2008) may inhibit 
star formation, perhaps driving the gas from the galaxy. 

The H I disks of galaxies, therefore, serve as a barome- 
ter of recent star formation activity and a record of prior 
processing. The mass, metallicity, velocity field, surface 
density profile, etc. reflect both the underlying dark mat- 
ter potential and also the star formation history of the 
galaxy. In the local universe, H I disks are mapped in the 
21cm line with radio telescopes. These data reveal the 
mass, surface density profiles, and kinematics of modern 
H I disks (e.g. Zwaan et al. 2005a; Walter et al. 2008). 
With current facilities, unfortunately, it is impossible to 
survey H I disks in 21cm emission at high redshift. Such 
analysis awaits the construction of facilities like the pro- 
posed Square Kilometer Array. 

In lieu of 21cm observations, one may observe H I gas 
through electronic transitions, e.g. the Lyman, Balmer, 
and Paschen series. Although these lines can be stud- 

1 Department of Astronomy and Astrophysics, UCO/Lick Ob- 
servatory; University of California, 1156 High Street, Santa Cruz, 
CA 95064; xavier@ucolick.org 

2 Department of Physics, and Center for Astrophysics and Space 
Sciences, University of California, San Diego, Oilman Dr., La Jolla; 
CA 92093-0424; awolfe@ucsd.edu 



ied in emission, they arise in ionized gas via recombi- 
nation processes. Furthermore, detectable fluxes require 
a strong (i.e. local) excitation/ionization source and one 
tends to map special, isolated regions of the galaxy. To 
study the bulk of an H I disk, one may instead probe 
the gas in absorption via the Lyman series 3 (Wolfe et al. 
1986). At the characteristic column densities of H I disks, 
the Lya transition is damped and astronomers refer to 
the observed profiles as damped Lya systems (DLAs; 
Wolfe et al. 2005). These DLA profiles are mainly re- 
vealed in the spectra of distant quasars, yet they also 
manifest in the spectra of GRB afterglows (e.g. Chen 
et al. 2005). Unfortunately, these intrinsically luminous 
sources cover only a small fraction of the sky such that 
one rarely intesects a given H I disk with multiple sight- 
lines (e.g. Ellison et al. 2007). Therefore, H I disks at 
high z must be studied statistically through the observa- 
tion of thousands of quasars across the sky. 

This experiment has been realized over the past few 
years as an unintended consequence of the Sloan Digital 
Sky Survey (SDSS) of high redshift quasars (Prochaska 
& Herbert-Fort 2004; Prochaska et al. 2005, ; here- 
after PHW05). In this Letter, we report on the results 
for a survey of the SDSS Data Release 5 (SDSS-DR5; 
Adelman-McCarthy et al. 2007). We place new con- 
straints on the projected H I column density distribution, 
total covering fraction, and integrated mass density of H I 
disks at z > 2. We search for evolution in these quantities 
from z = 2 to 4 and also compare the measurements with 
H I disks from the local universe (Zwaan et al. 2005b). 
These results offer new insight on the evolution of H I 
disks and their role in the processes of galaxy formation. 

3 Although 21cm absorption also traces H I gas, the optical 
depth is inversely proportional to temperature and one is primarily 
sensitive to cold, H I gas. 
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Fig. 1. — Left: The integrated frequency distribution f(Nm,X) of projected H I column densities for galaxies at z = 2 to 4. The 
overplotted curve represents the best-fit, double power-law which has a break column density = 10 2155 cm — 2 , a 'faint-end' exponent 
03 = —2.00 ± 0.05, and a high-end exponent 04 < —4.4 (95% c.L). (right) Cumulative distribution functions of f(Nm,X) for galaxies in 
a series of redshift intervals. For z > 2, the shape of /(;Y[u , X) is invariant and, remarkably, matches with the observed function for local 
H I disks (Zwaan et al. 2005b). The gray dotted line shows the results for z ~ 1 from Rao et al. (2006). 



Throughout the Letter, we adopt a ACDM cosmology 
with = 0.7, Q m = 0.3, and H = 72 km s^Mpc" 1 . 

2. RESULTS 

We have surveyed damped Lya systems at z > 2.2 
using the database of quasar spectroscopy from the 
SDSS-DR5. We have implemented the algorithms de- 
veloped in PHW05 to search for DLA candidates, mea- 
sure the survey path, eliminate strong BAL quasars, 
and to fit Voigt functions to candidate DLAs. Ta- 
bles and figures for the Lya fits are provided online 
(http://www.ucolick.org/~xavier/SDSS). The full sta- 
tistical sample now comprises 738 DLAs with Ajji > 
10 20 - 3 cm" 2 over a total redshift path Az — 3082.5, each 
located at velocity 5v > 3000 km s -1 from the back- 
ground quasar. 

In Figure la we present the Nm frequency distribution 
/(A^hi, X) of the full statistical sample. This measure de- 
scribes the projected column density distribution of H I 
gas in galaxies at high z per comoving absorption path- 
length dX . Similar to our previous results (PHW05), 
f(Nm,X) is well described by a power-law at low A^hi 
values f(Nm,X) ~ with a w —2, but transitions 



to a steeper function at Am ~ 10 215 cm 2 . This break 
in /(Ahi,A) is required to yield a finite H I mass den- 
sity, phi = (m p H /c) J N f(N m , X)dXdN . Following 
the formalism in PHW05, we fitted a double power-law 
to /(Ahi,A); the best-fit model is overplotted on the 
data in Figure la and tabulated in Table 1 for a series of 
redshift intervals. The table also lists the zeroth and first 
moments of f(Nm,X) which give the line density £(X) 
and phi values respectively. The former quantity rep- 
resents the covering fraction per dX for the integrated 
population of H I disks at a given epoch. The latter 
quantity is the comoving mass density of H I gas in high 
z galaxies. 

A principal result of our survey is that the shape of 
the H I distribution function does not evolve. Figure lb 
shows the cumulative f(Nm,X) distributions for a se- 
ries of redshift intervals from z = 2.2 to 5. We have 
performed a series of two-sided Kolmogorov-Smirnov test 
comparisons and find that the null hypothesis cannot be 
ruled out at greater than 90% confidence for any pair. 
Table 1 also reveals that the parameters of a double 
power-law fit to the data show no significant (> 2a) vari- 
ations with redshift. 
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a Mean absorption redshift of the DLA sample. 

b Line density of DLAs per absorption length dX. Also written as dN/dX in the literature. 

c Note that we recover the same estimate for phi whether we sum the discrete TVjji values or integrate the best fitting double power-law 
(see also PHW05). 



In Figure lb, we also present the cumulative f(Nm, X) 
function for H I disks in the local universe, as estimated 
from 21cm observations (Zwaan et al. 2005b). Remark- 
ably, the z ~ distribution function is a near perfect 
match to the high z universe. This is a stunning result 
noted first and independently by Zwaan et al. (2005b) 
and PHW05. Although the universe and the galaxies 
within it have evolved substantially over the w 10 Gyr in- 
terval from z = 3 to today, the combined distribution of 
H I surface densities is invariant. Stated differently, the 
population of galaxies at any epoch shows a self-similar 
projected H I surface density distribution. 

Although the shape of f(Nm,X) is invariant, its nor- 
malization decreases with time. This is revealed in Fig- 
ure 2 where we present £(X) and pui from z = 5 to 
2. Both the co-moving covering fraction and the mass 
density of H I disks decreases by 50% in this ~ 2 Gyr 
interval. This sharp decline in both £(X) and phi is 
a surprising and profound result. Before discussing its 
origin, we emphasize that the evolution must occur at 
all column densities of the DLAs such that the shape of 
f(Nm,X) remains invariant. Therefore, one should fo- 
cus on processes that affect the inner and outer regions 
of H I disks together. 

3. DISCUSSION 

The results of the previous section have far reaching 
implications for the nature and role of H I disks and for 
the processes of galaxy formation. Let us begin with 
the invariance in the shape of f(Nm,X). In terms of 
statistical power, the shape of f(Nm,X) is dominated 
by systems with low JVhi values and the primary result 
is that the 'faint-end slope' of f(Nm,X) is invariant. In 
the local universe, low A^hi sightlines correspond to the 
outer regions of H I disks. We draw the inference that 
galaxies have self-similar surface density profiles in the 
outer disk at all cosmic times. 

In the most straightforward, analytic models of galaxy 
formation (e.g. Mo et al. 1998), f(Nm,X) is determined 
by the radial H I surface density profiles, which in turn 
are set by the total mass of the system, the angular mo- 
mentum distribution of the galaxy, the gas mass frac- 
tion, etc. This simple picture is modified by spiral den- 
sity waves, warps, galaxy mergers, the detailed nature 
of ISM clumping, molecular cloud formation, and feed- 
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Fig. 2. — Upper: Comoving H I mass density phi of galaxies 
at z > 2 assuming a ACDM cosmology. The H I mass density 
is observed to decline by 50% from 2 = 4 to 2.2, an interval 
spanning less than 2 Gyr. The red band shows the estimate of Phi 
at z ~ from 21cm surveys of local H I disks (Zwaan et al. 2005a). 
Lower: The line density of DLAs per comoving absorption length 
dX. This quantity can be visualized as the integrated covering 
fraction per comoving pathlength for H I disks. Folowing the mass 
density phi 1 th e covering fraction decreases by 50% from z sa 4 to 2 
where it reaches the present-day value (red band) as estimated from 
21cm observations (Ryan-Weber et al. 2003; Zwaan et al. 2005a). 
Taken together, the results argue that H I disks have not evolved 
significantly over the past ?s 10 Gyr. 

back from supernovae and/or AGN activity. All of these 
processes are expected to vary with time, especially the 
characteristic mass of galaxies. The results presented in 
Figure lb suggest that the outer regions of H I disks are 
not especially sensitive to these processes nor to the un- 
derlying dark matter halo mass. We note that this is 
actually a prediction of viscous models of galactic disk 
formation (Lin & Pringle 1987; Olivier et al. 1991). We 
await explorations of this topic witin the context of cos- 
mological simulations of galaxy formation (e.g. Nagaminc 
et al. 2004; Razoumov et al. 2006; Pontzen et al. 2008). 
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Secular and feedback processes may be expected to 
have greater effect on the gas toward the inner re- 
gions, i.e. at the highest surface densities. We have also 
searched for variations in f(Nm,X) at large Am, but 
identify none. At the 95% c.L, all of the redshift inter- 
vals have /(Am, X) distributions consistent with a break 
column density of Nd = 10 21 ' 6 cm~ 2 . Furthermore, re- 
stricting the frequency distributions to Am > 10 21 cm~ 2 , 
all give satisfactory KS-test probabilities. We conclude 
there is no evolution at these column densities, but cau- 
tion that the full sample includes only 105 DLAs. The 
data also reveal, for the first time, that /(Am, A) is 
steeper than a = —3 beyond the break. This cannot be 
attributed to projection alone (which predictes a = —3; 
see Wolfe & Chen 2006), but we associate the steeper 
drop to the conversion of atomic gas to molecules (Schaye 
2001; Zwaan & Prochaska 2006). 

Now consider the sharp decrease in the total comoving 
covering fraction and H I mass density from z = 4 to 
2 (Figure 2, Table 1). One's initial reaction may be to 
interpret this decline in terms of active star-formation, 
i.e. the conversion of the H I gas in DLAs to stars via 
in situ star formation. This interpretation is problem- 
atic for several reasons. First, one expects star forma- 
tion to mainly influence gas at high H I surface densi- 
ties; this is revealed, in part, by the form of the Schmidt 
law, Ssfr = AxE 1 ' 4 , where Esfr is the SFR per unit 
area and X is mass surface density (Schmidt 1959; Ken- 
nicutt 1998). But the invariant shape of /(Am, A) sug- 
gests that in situ star formation in DLAs is unlikely to 
consume gas according to the Schmidt law. This is il- 
lustrated in Fig. 3a, which shows how an initial single 
power-law approximation for /(Am, A) steepens with 
time if stars form according to the Schmidt law (Lanzctta 
et al. 1995). We find that the absence of changes in the 
shape of /(Am, A) implies that the star-formation effi- 
ciency is less than 1/10 that in local galaxies. 

Second, Wolfe & Chen (2006) used the infrequent 
detection of extended, low surface-brightness galaxies 
in the Hubble Ultra Deep field to set an upper limit 
on the comoving SFR density of SFR density, p\ < 
10~ 2 - 7 M Q yr _1 Mpc -3 . We used this limit to set an up- 
per limit on the decrease in pm due to in situ star forma- 
tion. We fitted a conservative expression for p*(z) to be 
consistent with (a) this upper limit at z ~ 3 and (b) ob- 
servations of star-forming galaxies in the redshift interval 
z=[0,8] (e.g. Bouwens et al. 2008). We then computed 
the decrease in pm(z) by integrating p*(z) from z—6 to 
z. The resulting dip in pm is shown as the dotted curve 
in Fig. 3b. Clearly the decrease in pm predicted by in 
situ star formation is too small to account for the fac- 
tor of two decrease observed. Specifically, the total mass 
density of stars formed in-situ in H I disks is less than a 
few percent of the observed pm at z = 4. 

Third, to markedly change the covering fraction of the 
H I disks, in situ star-formation would have to lower 
« 50% of the H I disks below the DLA criterion. At 
the DLA threshold, star formation is likely very weak 
(if not absent) and should not affect this gas. We con- 
clude, therefore, that in situ star formation on its own 
is insufficient to explain the sharp decline in £(X) and 
Phi at high z. By similar arguments, one rules out an 
interpertation where the majority of H I gas is simply 
conveted into molecular gas. Furthermore, observational 



biases (e.g. dust obscuration, gravitational lensing) are 
most important at high Am values and can be ruled out 
as dominant factors. Finally, the extragalactic UV back- 
ground (EUVB) will modify the distribution of neutral 
gas in the outer parts of H I disks but the DLA threshold 
is sufficiently large that variations in the EUVB should 
play only a minor role for £(X) and a negligble one for 
Phi (e.g. Viegas 1995). 

On the other hand it is possible that the gas in DLAs 
fuels star formation (Wolfe & Chen 2006; Wolfe et al. 
2008). In that case the Schmidt law would not apply to 
DLA gas. Rather, secular processes could drive the gas 
to the center where it would be consumed by star for- 
mation in compact star- forming regions. To compute 
the decrease in pm(z) we fitted analytic functions to 
the intrinsic values of p*(z) derived by Bouwens et al. 
(2008) and computed pui(z) by integrating the under the 
p*(z) curve. The results are well-matched to the data at 
z=[2.5,4.5] (Fig. 3b). Although we obtained this fit by 
arbitrarily increasing the "initial" pm(z) at z—6 by 10 
% above our highest data point at z=4.0, the decrease in 
pni(z) predicted by this model for p*(z) provides a good 
match to the data for a wide range of initial values of 
Phi- However, the model breaks down at z < 2.4, where 
the predicted pni(z) falls below the constant value set by 
Pni(z) at £=[0,2.4]. But the levelling off of pm at z < 
2.4, could be explained by delayed infall of gas from the 
IGM at a rate that balances gas consumption by star 
formation, a phenomenon observed in some numerical 
simulations of galaxy formation (Keres et al. 2005). Of 
course, a crucial challenge to this idea is how to transport 
gas from the extended regions comprising DLAs to the 
center, i.e., by a factor of 15 or more in radius (cf. Wolfe 
et al. 2008), without affecting the shape of /(Am, A), 
but by reducing its normalization by a factor of two. 

The other extrcmum is that the sharp decline in cov- 
ering fraction and mass density of H I disks results from 
'violent' feedback processes. By violent, we envision pro- 
cesses that altogether remove the H I gas from a galaxy. 
These may include AGN activity, galactic-scale winds, 
tidal effects, and ram-pressure stripping. To match the 
observed evolution in £(X) and pm, one would require 
that approximately half of the galaxies exhibiting H I 
disks at z « 4 have lost their gas by z = 2. This implies a 
dramatic evolution in the fraction of H I-rich to H I poor 
galaxies in the 2Gyr interval centered at z « 3. The 
resulting galaxies, if unable to accrete new H I gas for 
subsequent star formation, would passively evolve into 
'red and dead' galaxies. We speculate, therefore, that 
z ~ 3 marks the formation epoch for the formation of 
stars in the ancestors of modern, early-type galaxies. 

The data in Figure 2 provides another surprising re- 
sult. The red bands in the figure show z « estimates of 
£(X) and pm, as inferred from 21cm observations (Ryan- 
Weber et al. 2003; Zwaan et al. 2005a). We find that p m 
at z = 2.2 matches the present-day value. A brazen, but 
reasonable, assertion by interpolation is that pm has not 
evolved over the past lOGyr of our universe. If con- 
firmed, this result has several important implications. 
First, if processes destroy H I disks below z — 2, these 
must be matched by the formation of new galaxies. Be- 
cause the assembly of dark matter halos of galactic-scale 
masses (< 10 12 M Q ) is expected to be nearly complete at 
z <~ 2, we contend that the destruction of H I disks is also 
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Fig. 3. — Models for the evolution of f{Nm,X) and pjjl at high z. The blue solid curves in the left panel show the evolution in /(-/Vhi, X) 
from z = 4 to 2 preidcted by assuming star formation occurs in situ according to the Schmidt law (Hsfr °c S^f) in a closed box model. 
This model predicts too much evolution in f(Nm,X). The smooth curves in the right-hand panel show the evoultion in phi predicted by 
the simple assumption that the H I consumption matches the mass of stars formed from z = 5 to 2 by integrating assumed star-formation 
histories. The green dotted curve shows the results for them situ star-formation history of DLAs consistent with the upper limits on p*(z) 
set for the Hubble Ultra Deep Field (Wolfe &; Chen (2006)). The green dashed curve shows the results for the star formation history of 
UV-selected galaxies (Bouwens et al. (2008)). This curve is arbitrarily normalized to the observed phi value at z = 4. 



nearly complete. Second, the accretion of gas into exist- 
ing H I disks must be balanced by the the consumption 
of that gas into stars and/or its removal from the galaxy 
by feedback processes. This suggests that H I disks play 
a special, but subserviant role in the formation of stars. 
In essence, H I disks represent a bias level of gas that 
is a necessary but insufficient condition for star forma- 
tion. In this scenario, the global star formation rate at 
any given epoch is driven predominantly by the accretion 
rate of fresh material onto existing H I disks, an infer- 
ence also drawn from cosmological simulations of galaxy 
formation (Keres et al. 2005). It further suggests that 
H I disks have been critically unstable to star formation 
over the past 10 billion years and possibly all of cosmic 
time. 

Before concluding, we comment briefly on current esti- 
mates of f(Nm, X) and pm at z ~ 1 where Mg II absorp- 
tion has been used to identify DLA candidates (Rao & 
Turnshek 2000; Rao et al. 2006). In Figure lb, we present 
the cumulative /(-/Vjji, A) distribution for DLAs derived 
in this fashion. The z ~ 1 sample exhibits a higher inci- 



dence of DLAs with AT H j > 10 21 cm~ 2 than at z > 2 or 
z ~ 0. If taken at face value, the results indicate that H I 
disks at z ~ 1 have a higher cross-section at column den- 
sities N m m 10 21 - 6 cm~ 2 than N m « 10 21 ' 3 cm~ 2 . It is 
unlikely that there is a single galaxy in the universe with 
this characteristic, much less the integrated population 
of H I disks. As such, we conclude that the z ~ 1 results 
suffer from a statistical fluke or an unfortunate observa- 
tional bias. If we demand that /(Ajji, X) at z ~ 1 follow 
the same shape as at z > 2 and z ~ 0, then we estimate 
a phi value consistent with no evolution over the past 
lOGyr. 

The observations presented here will be supplemented 
by future data releases of the SDSS and next generation 
surveys. The key open empirical issues include: (i) is 
there even a mild evolution in the break column den- 
sity with redshift? (ii) what is the functional form of 
f{Nm.iX) beyond the break? (iii) what are the values 
of these quantities at z m 6? Of greater interest will be 
to compare the observational constrainst presented here 
against theoretical models for the build-up and evolution 
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of H I disks over cosmic time. 
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